Universal shape ratios for polymers grafted at a flat 

surface 



Maria Serena Causo 
INFM-NEST and Scuola Normale Superiore 
Piazza dei Cavalieri, 7 
1-56100 Pisa, ITALY 
Internet: causoSsns . it 



February 1, 2008 



Abstract 

We consider dilute non-adsorbed polymers grafted at an impenetrable sur- 
face and compute several quantities which characterize the polymer shape: the 
asphericity and the ratios of the eigenvalues of the radius-of-gyration tensor. 
The results are only slightly different from those obtained for polymers in the 
bulk, showing that the surface has little influence on the polymer shape. 



1 Introduction 



It has been known since 1934 that polymers are not istantaneously spherical in 
shape 1^. Since then a lot of work has been devoted to quantify the average de- 
viation from spherical symmetric shape in the limit of large degree of polymerization, 
since this deviation plays an important role in the hydrodynamic behaviour of di- 
lute polymer solutions 0. There are several observables which quantify the shape of 
polymers. One considers the (squared) radius-of-gyration tensor and its eigenvalues 
qi, . . . ,qd, which give the square lengths of the polymer along the principal axes of 
inertia [^. The ratios {qi) / {qj) or, alternatively, (qi/qj) are universal and may be 
used to give a universal characterization of the polymer shape. Alternatively, one 
may use the shape asphericity that allows to distinguish between spherical and rod- 
like configurations. The equivalent ellipsoid has been experimentally determined in a 
small-angle x-ray scattering experiment for globular proteins 

The shape of Gaussian polymers has been extensively studied. The shape as- 
phericity, which in its two variants is proportional to the istantaneous variance of 
the average of eigenvalues of the radius-of-gyration tensor, has been exactly com- 
puted 0,11 . Instead, the ratios of the eigenvalues q^ are not exactly known. There 
are perturbative estimates at second order in 1/d 0, which are, however, not very 
predictive in the physically relevant case d = 3, and several estimates from numerical 
simulations [0-0. The numerical estimates are in good agreement with each other. 

In the case of polymers with excluded-volume interactions, asphericity has been 



computed to second order in e = 4 — ci in Ref. |]14[ and by numerical simulations 
in Refs. [^, while the ratios of the eigenvalues of the radius-of-gyration ten- 

sor have been computed by Monte Carlo simulations in Refs. p], |18|-Q . Here 
we consider the case of polymers with excluded-volume interactions grafted at an 
impenetrable interface at the ordinary transition p3|-p5[]. Polymers are modeled as 
A^-step self- avoiding walks on a simple cubic lattice. We perform an extensive Monte 
Carlo simulation using a variant of the pivot algorithm PB|-P5|, the cut-and-permute 
algorithm, which has been shown to be optimal for non-adsorbed polymers grafted at 
a surface The efficiency of the algorithm allows us to simulate long walks with 
N < 32000 and to study how the presence of the surface changes the shape of the 
polymer. 

We find that the main effect of the presence of the surface is to increase the 
length of the longest axis of the inertia ellipsoid with respect to the other two, whose 
relative length remains, instead, of the same order. Nevertheless, being this increment 
in asphericity only of a few percents, we expect that it can only be detected in high- 
resolution experiments. The paper is organized as follows. In Sec. 2 we define the 
observables that will be computed. In Sec. 3 we present our Monte Carlo results, 
paying due attention to the corrections to scaling. In Sec. 4 we compare our results 
with those for polymers in the bulk. 
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2 Definitions 



We model the polymer as an iV-step self-avoiding walk on a simple cubic lattice, 
each step representing a single monomer. The polymer has one of the two endpoints 
grafted at an impenetrable flat surface of equation z = 0, so that the accessible space 
is = {r = (ry, z) : ry G M^, z > 0}. Let us denote hy u = {uq, . . . , un} a generic 
polymer configuration, where Ui = {uj} , ujf , ujf) G Mi]_ is the position vector of the i-th 
monomer and = 0, and let us call Wiy the set of allowed configurations. 

One can introduce several invariant quantities which describe the shape of the 
polymer. In generic dimension d we define the radius-of-gyration or shape tensor 



whose trace concides with the squared radius of gyration 

RI = TtQ. (2) 
If we denote with {qi, q2, . . . , qd} the set of eigenvalues of Q, ordered in such a way 

Ir 
d 



that qi < q2 ^ ■ ■ ■ ^ qd, and with q = ^TrQ their average, we can define the shape 



asphericity of the polymer as 



where Q = Q — ^(TrQ) Id and Id is the identity matrix. Since the eigenvalues of Q 
correspond to the squared lengths of the equivalent ellipsoid axes and the eigenvectors 
represent the axes directions, Ad{u!) = for a spherically symmetric configuration — 
in this case all eigenvalues take the same value — while ^^(a;) = 1 if the polymer is 
completely elongated and all eigenvalues except one vanish. The mean asphericity is 
obtained by averaging Ad{uj) over the allowed polymer conformations 

A^ = -Y1 M^) , (4) 

where c^r is the number of allowed configurations and behaves as cjv ~ fi^N'''^"^ in 
the critical limit N —>■ oo. Here fi is the critical fugacity, which does not depend 
on the presence of the surface, and 71 is the susceptibility critical surface exponent 
7i = 0.679(2) (Ref. [§g), 71 ^ 0.680 (Ref. |HT|). Other invariant quantities are often 
considered in order to characterize the polymer shape. For instance, one may consider 

^ d (Tr(g^)) 

ci-l((TrQ)2) ' ^""^ 

or the relative magnitudes of the averages of the eigenvalues, {q2)/{(li), {iz) / {(li) ■ 
The advantage of studying these quantities is that they do not require computing an 
average of the ratio of highly fluctuating quantities. 
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In the limit of large degree of polymerization N ^ oo, these quantities reach 
universal values with corrections of order N~^, 

= A, + ^ + , (6) 

A^a = A.(l + ^ + ..^. (7) 
(^2"^) _ r2 A , ^2,1 , 



Lii 1 + 



where A is the leading correction-to-scaling exponent. On the other hand, the eigen- 
values of the gyration tensor scale as 

(gf ) = fe.iV^^ (1 + ^ + ■ ■ ■) . = 1,2,3. (10) 

In the following section we present a Monte Carlo study of the universal behaviour 
of these quantities. 



3 Monte Carlo study of the shape asphericity 

We computed the eigenvalues of the shape tensor and the asphericities by Monte 
Carlo simulations. The system was modeled as a self-avoiding walk grafted at the 
origin of the half cubic lattice Zi. In order to sample the phase space W^, we used a 



nonlocal Monte Carlo algorithm, the cut-and-permute algorithm |^9[. It is an optimal 
algorithm for studying polymers grafted at an impenetrable surface at the ordinary 
transition, since autocorrelation times in CPU units for global observables are simply 
proportional to the degree of polymerization A^. The algorithm has, therefore, the 
same optimal behaviour as the pivot algorithm for polymers in the bulk. 



The system in the bulk is known to have strong corrections to scaling p2| 



smce 



the leading correction-to-scaling exponent is Abuik ~ 0.5. It is expected that polymers 
in the presence of a surface show analogous corrections, and, indeed, our estimates of 
the leading correction-to-scaling exponent are close to 0.5. 

Since the presence of large corrections to scaling introduces systematic errors in 
the estimates for finite A^, we studied the system for large A^, 1000 < N < 32000. At 
first, we studied the behaviour of the eigenvalues of the shape tensor and of its trace 
i?^. The results are displayed in Table 0. These quantities are expected to increase 
as A^^^, being v the same critical exponent as in the bulk. In order to check this 
prediction, we performed fits of the form 



(g,)jv = hN''' ^ = 1,2,3, 

{Rl) = bN'\ (11) 
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where b, bi, and u are free parameters. In Eq. (|TT]) corrections to scaling are neglected. 
We expect, therefore, a systematic error in our estimates. In order to take into account 
the effect of corrections to scaling, we have performed different fits, considering only 
data with N > A^min- The effective exponents which are obtained decrease with 
increasing A^min for all values of A^min considered, so that we can only give an upper 
bound for the exponent u. The estimates which are less affected by corrections to 
scaling appear to be those corresponding to the largest eigenvalue. Therefore, taking 
as an upper bound for z/ the estimate from the fit of (q^) with A^min = 16000, we obtain 
u < 0.58774(35), which agrees with the best estimates obtained so far for polymers 
in the bulk, which are z/ = 0.58758(7) (Ref. H), u = 0.5877(6) (Ref. |5[). 



Let us consider now the ratios of eigenvalues. For such invariant quantities we 
performed the following three-parameter fits 

^"■^"-LL + Df^Y k = l,2, (12) 



where i, D, and A are free parameters and A^o = 16000 is a fixed parameter intro- 
duced in order to reduce the cross-correlations between D and A in the correction- 
to-scaling term. Again, since we are neglecting higher-order corrections to scaling, 
the estimates of the universal amplitudes are affected by systematic errors. In order 
to take into account this effect, we performed fits in which only data with A^ > A^min 
were considered and we studied how our estimates approach the asymptotic value as 
A^min increases. The results of the fits are shown in Table ^. In this case the esti- 
mates from the three-parameter fits appear to be very stable, so that we can choose 
as best estimates of the universal ratios the ones for A^min = 1000: = 2.9272(28), 
L|i = 14.6293(66). 

Nevertheless, one can notice that, while the correction-to-scaling exponent deter- 
mined from the ratio {q3)N/{(li)N is A = 0.450(46), which is in agreement with the 
estimates for the leading correction-to-scaling exponent for polymers in the bulk , 
the estimate of A obtained in fits of {q2)N/{qi)N is considerably smaller. This can 
be due to the presence of corrections with opposite sign and could affect our final 
estimate. For this reason, we also performed fits in which the leading correction-to- 
scaling exponent was fixed at A = 0.5, and higher-order corrections were neglected. 
The estimates obtained are shown in the last three rows of Table ^. Estimates for i 
are stable for A'min = 4000, while those for L| ^ are already stable for A^min = 1000. 
In all cases the P^r degree of freedom is less than 1. Therefore, we can take as 
final estimates 

= 2.94931(99), (13) 
Ll^ = 14.6449(60). (14) 



One can notice that the estimate in Eq. (13) agrees with the one obtained without 
fixing A, while 1 in Eq. (0) is not compatible with the previous estimate. Since we 
know a priori that A should be approximately 1/2, we are inclined to consider ( p!^ 
a more reliable estimate. 
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We performed the same kind of analysis for the asphericities, performing fits of 
the form 

< = A, + D(^^y, (15) 

= ^3+5(f)^ (16) 

The results are reported in Table |^. Fits are not very stable and the amplitude 
appearing in the correction to scaling is very small. For A^min = 1000 we estimated 
D ^ {47 ± 10) X 10"^ and D ^ {43 ±7) x 10"^. As before, we also performed a 
two-parameter fit in which the leading correction-to-scaling exponent was fixed at 
A = 0.5. Since all estimates for iVmin > 1000 are compatible within error bars, our 
final estimates are 

As = 0.44520(15), (17) 
As = 0.55949(13). (18) 

These estimates are in agreement with those obtained letting A be a free parameter. 



Comparison of the results with the case of poly- 
mers in the bulk 



14,18-22 



Let us compare our results with the ones for polymers in the bulk 
An early Monte Carlo simulation of polymers with length < 220, in which the 
presence of corrections to scaling was taken into account, gave A3 = 0.431(2). Since 
then, more computing power and more efficient algorithms became available and much 
longer polymers could be studied. In Ref. [|l^], self-avoiding walks for a single value 
of (A^ = 963) were simulated using the pivot algorithm, obtaining ^3 = 0.43333(8) 
and A3 = 0.54722(9) for the shape asphericities. From the data in Ref. for 
the shape factors, we can compute the eigenvalue ratios L| ^ = 14.1102(52) and 
L21 = 2.9631(12). Since only one value of A^ was considered, no analysis of the 
corrections to scaling could be done. Thus, we expect the errors to be grossly under- 
estimated. Indeed, our results show that, for A^ ^ 10^, the systematic error due to the 
scaling corrections is much larger than the quoted statistical error. In Refs. |]2ll,^| 
even longer walks were considered and data were extrapolated to infinite length with 
leading correction-to-scaling exponents A = 0.5 and A = 0.46. The estimates for 
the asphericities are: ^3 = 0.4302 and A3 = 0.5450 for A = 0.5; A3 = 0.4303 and 
A3 = 0.5452 for A = 0.46. From the asymptotic values for shape factors, we obtain 
for the eigenvalue ratios: L^^ = 13.89, L^^ = 2.934 for A = 0.5; and = 13.92 
= 2.939 for A = 0.46. No error bars on the extrapolated values were reported 
in Refs. [^^. The asphericity A3 has also been computed in an extensive Monte 
Carlo simulation in Ref. |18[ , where walks of length A^ > 40000 were considered. If 
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Table 1: Monte Carlo results for the eigenvalues and for the trace of the radius-of- 
gyration tensor Q and effective exponents computed using the results for walks of 
length > A^min = ^ ■ Data are fitted as in Eq. (|r^). "DF" is the number of degrees of 
freedom. 
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Table 2: Monte Carlo results for the ratios of eigenvalues and results of the corre- 
sponding fits for different values of A^min- Data were fitted to the scaling form (|12D. In 
the first set of fits A was taken as a free parameter, while in the second set A = 0.5. 
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Table 3: Monte Carlo results for asphericities and results of the corresponding fits for 
different iVmin. Data were fitted to the scaling form (|16]). In the first set of fits A was 
taken as a free parameter, while in the second set A = 0.5. 
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the data reported there are fitted assuming A = 0.5, one finds A3 = 0.43074(57), 
which agrees within error bars with the estimate for = 40000, ^3 = 0.4309(17), 



and with the above- reported estimates of Refs. [^,^. The same analysis cannot be 
performed for the eigenvalue ratios, since the estimates at different N are not avail- 
able. Nevertheless, using Table 2 of Ref. we obtain L| ^ ^ 14.05, ^ 2.963. 



Both values are slightly higher than the corresponding ones in Refs. pi|,p^. 

Comparing our estimates (P^, (|TBp, with those reported above, we find that, due 
to the presence of the surface, ^3 approximately increases by 3.5%, while A^, by 2.6%. 
The change in the two asphericity observables is, therefore, of the same order. 

Moreover, while ^ does not appreciably change, L3 increases by about 5%. 
This shows that the change in asphericity is mainly due to the fact that the length of 
the longest axis of the inertia ellipsoid increases with respect to the other two. This 
is a signal of the fact that the presence of the surface gives the polymer directionality 
along the normal direction to the surface, which we expect, on average, to coincide 
with the direction of the larger axis in the inertia ellipsoid. Nevertheless, the difference 
is only of a few percents. Therefore, we expect that it can be experimentally detected 
only in high-precision experiments. 
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